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Quantum fluctuations of an inflaton fleld, slow-rolling during inflation are coupled to metric 
fluctuations. In conventional four dimensional cosmology one can calculate the effect of scalar metric 
perturbations as slow-roll corrections to the evolution of a massless free field in de Sitter spacetime. 
This gives the well-known first-order corrections to the field perturbations after horizon-exit. If 
inflaton fluctuations on a four dimensional brane embedded in a flve dimensional bulk spacetime 
are studied to first-order in slow-roll then we recover the usual conserved curvature perturbation on 
super-horizon scales. But on small scales, at high energies, we find that the coupling to the bulk 
metric perturbations cannot be neglected, leading to a modified amplitude of vacuum oscillations 
... on small scales. This is a large effect which casts doubt on the reliability of the usual calculation of 

' inflaton fluctuations on the brane neglecting their gravitational coupling. 

o ■ 

O ■ PACS numbers: 04.50.+h, 98.80.Cq 

^ ■ I. INTRODUCTION 

^ ; 

CO ' Inflation is probably the simplest scenario for the origin of primordial fluctuations in our Universe 0]. Small scale 
vacuum fluctuations can be stretched to astrophysical scales by an period of accelerated expansion. Inflation provides 
a test of high-energy physics because the perturbations are generated from very short scales at high energies in the very 
^ , early universe. These perturbations carry signatures from high energy physics, which can be tested by astronomical 
I observations. 

The slow-roll approximation 2] is a useful tool to study the fluctuations generated during inflation. If we can 
neglect the coupling to metric perturbations and the effective mass of the field then the perturbations are described 
by the fluctuations of a free scalar fleld in de Sitter spacetime. This gives the familiar result that the power spectrum 
■ of scalar field perturbations at horizon-crossing is given by {H/2tt)'^. One can then calculate the comoving curvature 
perturbation which is conserved on super-horizon scales for adiabatic perturbations. 

However infiaton perturbations will be coupled to gravity (metric perturbations) at first-order in the slow-roll 
parameters. In four-dimcnsional general relativity it is known how to consistently include linear metric perturbations 
by working in terms of the gauge invariant combination of scalar field and curvature perturbations, the so-called 
Mukhanov-Sasaki variable, which obeys a simple wave equation 3]. Gravitational effects are negligible at small 
scales and high energies, where perturbations can be normalised to the usual Bunch-Davies vacuum state. On large 
scales (super-horizon scales) the comoving curvature perturbation is conserved allowing one to relate observations of 
temperature anisotropics in the cosmic microwave background to high energy vacuum fluctuations during inflation. 
K> Exact solutions are known for the special case of power-law inflation in general relativity which generalise the de 
j_] Sitter result and have been used to calculate first-order slow-roll corrections in more general inflation models 0|. 

In this paper, we develop a new way to derive slow-roll corrections based on a slow-roll expansion about de Sitter 
spacetime. In four-dimensional general relativity we show how to recover the usual flrst-order slow-roll corrections. 
Our method may be useful when one cannot derive an exact solution and the background spacetime is given as a 
perturbation about de Sitter spacetime. 

We then apply our method to inflation in the brane world model. New ideas in the string theory suggest that our 
observable universe is a 4-dimensional hypersurface, or brane, in a higher dimensional bulk spacetime The simplest 
example of this model is the Randall-Sundrum model where there is a brane embedded in a 5-dimensional anti-de 
Sitter (AdS) spacetime |^. An AdS spacetime has a characteristic curvature scale /i associated with the negative 
cosmological constant in the bulk. The spacetime shrinks exponentially away from the brane and this geometry 
effectively compactifies the 5-dimensional spacetime with the effective size On large length scales L > 

4-dimensional Einstein gravity is recovered, while on small scales, the gravity becomes 5-dimensional 0- In the early 
universe when the Hubble horizon is smaller than /i^^, we expect significant effects from higher dimensional bulk 
spacetime. Indeed, the Friedmann equation is modified from the conventional 4-dimensional theory for H^i 2> 1 0. 
This modification of Friedmann can provide a novel model for inflation ilTol ITU . 

In Ref.Q, the amplitude of the curvature perturbation is calculated by taking into account the modiflcation of 
the Friedmann equation. This work has been extended to include higher order corrections in slow-roll parameters 
[l^ and the formula has been widely used to confront this model with the observations 13]. But to derive these 
formulae for the spectrum of the primordial curvature perturbations the effect of coupling to five-dimensional gravity 
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has been neglected and in particular it is assumed that the power spectrum of inflaton perturbations at horizon 
crossing is given by (iJ/27r)^. This assumption is only valid to zeroth order in slow-roll parameters. At first order 
the inflaton perturbations will be coupled to metric perturbations. In the brane world, metric perturbations live 
in the 5-dimensional spacetime, and thus we must check if 5-dimensional effects change the result of conventional 
4-dimensional theory. Especially, at small scales/high energies, the 5-dimensional effects could be large. 

The first attempt to study the backreaction due to metric perturbations was made in Ref p^ . There, perturbations 
are solved perturbatively in slow-roll parameters. We should emphasize that this is the only possible way to perform 
the calculations analytically. If the background spacetime of the brane deviates from de Sitter spacetime, we cannot 
solve the bulk metric perturbations analytically. In contrast to four-dimensional general relativity, there are no 
other exact solutions known for the perturbation equations. Thus we must develop a new approach to calculate 
the effect of slow-roll corrections. In this paper, we extend earlier studies and investigate the backreaction due to 
higher-dimensional perturbations using a slow-roll expansion. 

The structure of the rest of the paper is as follows. In section II, we describe our new approach to derive first order 
slow-roll corrections in a conventional 4-dimensional cosmology. In section III, we review an inflation model in the 
Randall-Sundrum brane-world driven by an inflaton field on the brane. In section IV, we derive the equations that 
govern the coupled system of inflaton fluctuations on the brane and metric perturbations in the bulk. In section V, the 
first order corrections to the inflaton fluctuations on the brane are solved. In section VI, we discuss the implications 
of our result for the brane world inflation model. 



II. SLOW-ROLL EXPANSION OF SCALAR PERTURBATIONS IN 4D COSMOLOGY 



A. Background spacetime 



We consider an inflaton cj) whose potential energy density V{4') drives inflation. In the conventional 4-dimensional 
general relativity described by the metric 

ds^ = -dt^ + a{tf5ijdx^dx' , (1) 
the Friedmann equation and the equation of motion for the homogeneous field, cj), are given by 

+ 3ff^--^, (3) 

where H — a/a, K4 — SnGi and G4 is the 4D gravitational coupling constant. A dot indicates a derivative with 
respect to cosmic time, t. Slow-roll parameters are defined by 

Slow-roll inflation is described by small values of e and rj. 



B. Slow-roll corrections to inflaton fluctuations 



The inhomogeneous inflaton fluctuation, 5(j), is coupled to the metric perturbations. In the Longitudinal gauge, the 
perturbed metric is written as 

ds^ ^ -(1 + 2m)dt^ + a{tf{l + 2^)S,jdx'dx'. (5) 
The coupled equations for S(f>, ^' and <f> can be simplified by using Mukhanov-Sasaki variable defined by 

u^a(^S^-^^y (6) 

Expanding u by Fourier modes, the wave equation for u is given by 
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where z = {acj)) / H and t is a conformal time defined as 

dt 



a{t) 

In the case of the slow-roll inflation, the mass term in Mukhanov-Sasaki equation Q can be approximated as 

ig = ^(2 + 6e-3, + 0(,^.^)), (9) 
up to first order of the slow-roll parameters. Then Eq. Q can be expressed as 

(Puk . ( ,1 1 



dT2 



fc2-^(2 + 6e-3?7) Ufe = 0. (10) 



Usually, the appropriately normalized solution with the correct asymptotic behavior at small scales is obtained by 
solving Eq. l(Tn|) directly as 

Uk{T) = ^e»(-+i/2)-/2(-r)i/2i/(i)(-fcr), (11) 

where v = 3/2-|-2e — r] and H'iF' is the Hankel function of the first kind of order u. Here we assumed the Bunch-Davies 
vacuum state where perturbations stay in Minkowski vacuum at small scales. Equation is an exact solution of 
the perturbation equation H10|l only if the slow-roll parameters e and ry are constant. However their variation in a 
Hubble time is second-order and hence of higher-order in the slow-roll expansion. Thus we can take e and rj to be 
evaluated around the time of horizon-crossing. 

We are interested in the asymptotic form of the solution well outside the horizon. Taking the limit — fcr — > yields 
the asymptotic form of Uk] 

r(3/2)V2fc^ ' ^ ' 

Expanding the gamma function in Eq. H12|l , we get 

Uk ^ eK-i/2)V2|i + _ ^)(2 _ ^ _ ln2)|^(-fcr)-i-2.+^, (13) 

J v2fc 

where we have used the formula for the poly-Gamma function 

V^(3/2)^^^ = 2-7-21n2, (14) 

where 7 is an Euler number. 

The quantity that is related to observables today is the the power spectrum of the curvature perturbation given by 
From Eqs. @ and H10|l . it can be shown that, at large scale, the time dependences of Uk and z are the same, 

1 /2 

that is, is constant. Note that the constancy of TZ in the large scale limit does not depend on the slow-roll 

approximation, but holds for any adiabatic perturbation. Thus this comoving curvature perturbation can be related 
to the perturbation in the radiation density on large scales long after inflation has ended. 
For the model with a monotonous potential, the following relation holds: 



dH 



" H ^ kIH 

and conformal time can be evaluated up to the first order in slow-roll parameters as 



(16) 



'-^<' + * (17) 

Thus the power spectrum of the curvature perturbation is given by 

where C = — 2 + ln2-|-7 ~ —0.73. The terms proportional to the slow-roll parameters are called Stewart-Lyth 
correction 
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C. Perturbing about de Sitter spacetime 

In this subsection, we reproduce the usual slow-roU corrections in a perturbative approach which does not require 
any exact solution of the perturbation equation other than that in a de Sitter spacetime. This will be more suited to 
extension to the case of brane-world gravity. 

At zeroth order in slow-roll parameters, the spacetime is described by the de Sitter spacetime. Thus we can expand 
the spacetime from de Sitter spacetime. The scale factor is expanded as 

a(t) =a(")(t)+a(i)(t)+0(e2), (t) = exp(fft). (19) 

Accordingly, the Mukhanov-Sasaki variable is expanded as 

Ufc(r)=4°)(r)+z.«(r)+0(e2), (20) 

where 4"^ = a<5<?!>(0) and = a((5(/)'-^' — {4'IH)'^). Substituting Eq. if^ into Eq. (fTn|l . the zeroth order equation is 
given by 



Since we expect that the effects of the deviation from de-Sitter spacetime are insignificant at small scales, the form 

ih-Davies vacuum 

u^°\T)^A{-Tf'^Hf/^{~kr), (22) 



of u"^^ is determined by demanding a Bunch-Davies vacuum 



where A = {^I2)e™ and the phase 6 is fixed so that Uk(T) (l/v2fc)e . 
Next, we must solve u]^^ sourced by this zeroth order solution 

^^+(fc^-A)4^)-l(6e-3rf^0. (23) 

If we impose the boundary conditions (i) u'^^t) is negligible in the limit t — > —oo, and (ii) M^^^(r) does not diverge 
faster than (t) in the limit t ^ 0, then we find that the solution is given by, 

4'^ = Ci(-fcr)i/2j3/2(-fcT)+C2(-fcr)i/2i/W(_fcr), 
Ci ^ ^(6e~3rMf dr'i {h« (-fcr')}' , 

C2 = ^(6e-37i)Ar dT'^Hi%{-kT')J3/2{-kT'), (24) 

where Ji, is the Bessel function of the order i/. 

We take the limit —kr—^O and compare the asymptotic form with Eq. (|13() . Using the small arguments limit of 
the Bessel functions 



•h/2 



r(5/2)^ 

we can show that the zeroth order Mukhanov-Sasaki variable approaches to 



uf^(T) -4z^^v4(2fc)-i/2(-fcr)-i. (26) 

TT 

Next, we must evaluate the asymptotic form of the first order Mukhanov-Sasaki variable. Using Eq. (j2Hl, Ci is 
evaluated as 

C, ^ ,U2.-,„a'-^j^^ (27, 
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We should be careful in evaluating the asymptotic behavior of C2 because sub-leading terms are comparable to the 
contribution from Ci. Using 

/ 2 / sin X \ / 2 / cos x\ 
J3/2{x)^\ — cosx , J_3/2(a;) = -W — smxH , (28) 

\ TTX \ X ) \ TTX \ X / 

the integral for C2 in Eq. H24|) can be evaluated as 

dr'^H^,%{~kT')J3/2{~kT') _|ici(-2fcr) + ^, (29) 
where Ci is the integrated cosine function defined as 

Ciix) = - r —dt. (30) 

Jx t 

For small — fcr, the integrated cosine function can be expressed as 

Ci{~2kT) ^ 7 + In 2 + In(-fcr). (31) 
Therefore, the asymptotic form of C2 is given by 

C2 -> (2e - r])A (^7 + In 2 + ln(-A:T) + . (32) 

Then we obtain the asymptotic form of Uk for — fcr — > up to the first order in slow-roll parameters 

Ukir) (-^)e'''|l + (2e-r?)(2-7-ln2)}|l- (2e-7?)ln(-fcT)}-i=(-fcr)-\ (33) 

where we have used the fact that A = (Y^/2)e'^. This should be compared with Eq. (|13|l . There appears a logarithmic 
term which diverges for —fcr 0. However, if we can renormalize this divergence by rewriting the logarithmic term 
as 

1- (2e-7;)ln(-fcT) ~ (-fcT)-2^+''. (34) 

we see that Eq. (|33ll is consistent with Eq. (|13() . 

Indeed, the logarithmic divergence for —fcr ^ in Uk does not show up in the spectrum of the curvature pertur- 
bation. In order to see this, we expand the curvature perturbation as 

V^' = {V^'Y"^ + {V^''}'^''> + Oie^). (35) 

On the other hand, by the definition of the curvature perturbation IjlSI) . we can write the spectrum of curvature 
perturbation up to the first order in slow-roll parameters as 



(36) 



where we also expanded z = {a4)) / H as 

^ = z(0)+z(i)+O(e2). (37) 

Since there is a difficulty to define the curvature perturbation in de Sitter spacetime, we concentrate on the ratio 
between the zeroth order and the first order of the curvature perturbation. By comparing Eq. H35|) to (|36|l . the ratio 
is given by 

{pV2}(0) „(o) ,(0)- ^ ^ 
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In order to evaluate Eq. H38|l . we must obtain z^^^ and z^^\ that is, we must solve Eq. ® perturbatively. Substituting 
Eq. H37|) into Eq. ©, the equation for z at zeroth order is given by 

If we consider only the growing mode, the zeroth order solution for z^^-* can be obtained as 

z(°)=Bt-\ (40) 
where B is an integration constant. This zeroth order solution gives a source term in the equation for z at first order; 



d^i) _2_ (66 -3r;) 



(41) 

The growing mode solution for the first order z*^^) is given by 

z^^) = -(2e - ?7)Bt"1 In(-fcT) + BDt~^, (42) 
where D is another integration constant. Then we get 

^ = -(2e-77)ln(-A;T) + i?. (43) 
This logarithmic divergence term exactly cancels the logarithmic divergence term in Uk] 

^ = (2_^_ln2-ln(-fcr))(2e-77). (44) 

From Eqs. (021 and iBH) we obtain 

{pV2}(0) 



-C(2e - ?/) - i?, (45) 



where C = —2 + ln2 + 7 ~ —0.73 is again a numerical constant. We cannot determine D in this approach, which 
comes from the difficulty to define curvature perturbation in pure de Sitter spacetime. However, we can still fix D as 
follows. Neglecting the logarithmic term, which is canceled by the contribution from Uk, the solution for z is written 
as 

z ^ B{1 + D)t-'^ . (46) 

This must be compared with the definition of z 

z^^M^jA^il^,)r~\ (47) 

where the solution for a up to the first order was used. Then we can identify B = — and D = e. Then Eq. I|45() 

agrees with the Stewart-Lyth correction given by Eq. (fl^ . 

III. SLOW-ROLL INFLATION IN RANDALL-SUNDRUM BRANE WORLD 

In this section, we apply our perturbative approach to the brane-world model. We consider the simplest version of 
brane-world inflation model based on the Randall-Sundrum model. We will consider a single brane embedded in a 
5-dimensional AdS spacetime. We assume that the inflaton is confined to the brane while gravity can propagate in 
the whole 5-dimensional spacetime 9] . 

The 5-dimensional metric describing this model is given by |^ 

ds^ = dy^ - N{y,tfdt^ + A{y,tf6^jdx''dx^, (48) 
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where 

A{y,t) = a{t) 



cosh fiy — ( 1 + — ) sinh 



N{y,t) = cosh^y- ^1 - ■^(2 + 3u;)^ sinh^y, (49) 

p^i^^ + Vicb), P^^<p'-V{q^), (50) 

and w = P/p. The brane is located at y = and the inflaton is confined to this hypersurface. On the brane, the 
Friedmann equation and the equation of motion for the scalar field are given by 

(51) 

* + 3H* = -^, (52) 

where = SttGs and is 5D gravitational coupling. We can define slow-roll parameters in the same way 

as the conventional cosmology, Eq. |0J. 

Unfortunately, the background metric H48|l is not in general a separable function with respect to y and t. Thus we 
cannot solve the metric perturbations analytically. In order to solve for the y-dependence of the bulk gravitons and 
to study the time-dependence of the perturbations on the brane, we will expand about the special case of a de Sitter 
spacetime on the brane. This corresponds to the background solution to zeroth order in a slow-roll expansion. For a 
de Sitter brane, AdS bulk gives a separable form for the bulk metric [Tsj : 

ds^ = dy^ + N^{y) [-dt^ + a^{t)5^jdx'dx'] , (53) 

where 

a{t) = (54) 

Niy) = -sinh/x(yh- I2/I), ■ (55) 
A* 

and y = ztj/jj are Cauchy horizons [l5j |. with 

y, = -coth-|^yi+^^j j • (^6) 

It is often useful to work in terms of the conformal bulk-coordinate z ~ j dy/N{y): 

z = sgn(y)i/-i In [coth \p{yx, - \y\)\ . (57) 
The Cauchy horizon is now at \z\ =00, and the brane is located sd z — ±Zb, with 

Zb = ^sinh-i-. (58) 
H p 

The line element, Eq. (|53|) . becomes 

ds^ = N^{z) [-dt^ + dz^ + c'^"'dx^] , (59) 

where 

N{z) = . ,. . (60) 

^smh(ii|z|) 
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IV. EQUATIONS FOR BULK METRIC PERTURBATIONS AND INFLATON PERTURBATIONS ON 

THE BRANE 

In this section, we derive the basic equations for the coupled Mukhanov-Sasaki variable on the brane and bulk 
metric perturbations following Ref . . 

A. Master variable for perturbations in AdS bulk 

In the background spacetime given by Eq. H59|l bulk metric perturbations can be solved using the master variable 
[id llTjl . The perturbed metric is given by 



ds^ = N{zf 



(1 + 2Ayy)dz^ + 2Aydtdz - (1 + 2A)dt'^ + 0^(1 + 2TZ)6ijdx'dx^ 



(61) 



In the special case of a de Sitter brane in the AdS bulk, the metric variables are written by the master variable Q as 



A 



An, 



' ( 2n" - 3—n' + n- fi^N^n 

6 \ N 

n" - 3 — Q' + 2n + fi^N^n 

6 V N 



R = 



6 



From the perturbed 5-dimensional Einstein equation, we can derive the equation for f2 



N' 



n - 3Hn - [n" - s—n' ] + —n - fi^N^n = o 



fc2 



(62) 
(63) 
(64) 
(65) 

(66) 



Solutions of the master equation can be separated into eigenmodes of the time-dependent equation on the brane and 
bulk mode equation: 



n{t,y;x) = / d^kdmam{t)um{z)t 



where 



Ctrri — 3Ha 
N' 



m 



3— uj„ + [i^N'^Um + vr?u„i = . 



(67) 



(68) 



Note that the Hubble damping term SHam has the "wrong sign", i.e., this is not the standard wave equation for a 
scalar field in four-dimensions. 

If we write am = a?^m and work in terms of the conformal time r = —l/{aH), the time-dependent part of the 
wave equation (|67|l can be rewritten as 



dr^ 



(^m = . 



This is the same form of the time-dependent mode equation commonly given for a massive scalar field in de Sitter 
spacetime. The general solution is given by 



ly = 



i?2 ' 



(69) 



where B^, is a linear combination of Bessel functions of order v. The solutions oscillate at early-times/small-scales 
for all m, with an approximately constant amplitude while they remain within the de Sitter event horizon (fc ^ aH). 
'Heavy modes', with m > ^H, continue to oscillate as they are stretched to super-horizon scales, but their amplitude 
rapidly decays away, oc a~^. But for 'light modes' with m < the perturbations become over-damped at 
late-times/large-scales {k <C aH), and decay more slowly: |u^J cx 0^"^'^. 
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B. Mukhanov-Sasaki equation on the brane 

Now we introduce a scalar field fluctuation on the brane. We expand the scalar field perturbation in terms of 
slow-roll parameters; 

50 = ,50(")+<50(i) + ... (70) 
The 0-th order of the scalar field fluctuation obeys the following equation of motion, 

+ 3H6<}>^°^ + ^Scp^"^ = 0. (71) 

The metric perturbations are generated by the 0-th order fluctuation of the scalar field through the induced Einstein 
equations on the brane [T^ . 

(72) 
(73) 



3iJ* - 


3iJ2$ + 






50o + 1/'<50W)H 


4*.. 






-ii = 




^^"^ - ^SqE, 




* - 3i?* + i/<i> + 3i72$ 


IP. 
3 


+ $) = 




5'(/)o-y'(50W)H 








+ $) = 


= kIStte, 







(74) 
(75) 
where 

A{y = 0,t) = <t>{t), R{y = 0,t) = ^{t) (76) 

45PE = (77) 
'^ISqe = ^^(Ji-Hri), (78) 



<5^E = —[n-Hn + ] , (79) 



and 

N' 

K4,cff = -Hi — 



(80) 



y=o 



The contributions SpE,SqE and (Jvtb come from the projected 5D Weyl tensor and these describe the effect of the 
bulk gravitational perturbations _1£] . The metric fluctuations in turn affect the dynamics of the first order scalar field 
perturbation 

/ 2 

S(f>^^^ + 3HS(j)^^^ + ^S(f>^^'> = -V"d(l>^"'> - 30* + 2y'$. (81) 

In order to evaluate the effect from metric perturbations, it is useful to use Mukhanov-Sasaki variable Q as in the 
conventional cosmology; 

Q^5cf>-^^. (82) 

In terms of slow-roll expansion, we have Q(o) ^ j0(o) ^^^^ q(i) ^ ^^(i) _ {^/H)^. The n using the induced Einstein 
equations, Eas. H72|l . H74|l and (|75(l . we can derive the equation for Q*-^^; 

Q(i) + 3i/Q(i) + ^Q(i) = -y"Q(") - 6ijQ(") -t- J, (83) 



10 



where 



J = 



H 6 



(84) 



The equation is the same as the standard 4-dimensional cosmology except for the term J, which describes the 
corrections from the 5-dimensional bulk perturbations. Because J contains the 5-dimensional quantity we must 
solve the bulk equation for SI to evaluate the effects. 

C. Boundary condition for SI 

In order to solve fi, we must specify the boundary condition for fl. we rewrite the expressions of $ and 4', Ea. H62() 
and inni, as ^ 



* = 



1 AT~3 



6 



6 



where 



N' ■ 

N' ■■ ■ 

N 



N' 

T = n' n. 

N 



(85) 
(86) 

(87) 



Substituting these expressions into the induced Einstein equations (|72|I - H75() . we obtain the equations written only by 
T and (50(0): 

- - k'a-^T = K2a(#0(o) + F'(</))50(°)) , 
= K^a{^5(i)^^) -V'{(j))6(t)^^'^). 



(89) 
(90) 



These equations can be thought as the boundary conditions for Q,. Combining the junction conditions, Eas. H88|) - (|90|l . 
we get an evolution equation for 



(91) 



This is consistent with the equation for scalar field equation Ea. H71|l . 



V. PERTURBATIVE SOLUTIONS 



We must solve the coupled equations Eqs. (|66|l for SI and Eq. H83|) for Q. Introducing dimensionless quantities 

Q(t) = Raify-^uij), fi(z, t) = K^^H'^uj{z, r), (92) 
the coupled equations are written as 

1 



,00/ 4 , o \ „ cosh Hz 

k^T^ I w + -cO + /c^tj I = uj" + 3 



sinh Hz sinh^ i7z 



= aH'^u, = iuj' 



cosh Hz 
sinh Hz 



k^u- —{2 + 6e-3r])u ^ Ju, Ju^-P^k^r^ 



UJ H LU + K UJ 



(93) 
(94) 
(95) 
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where a dot denotes a derivative with respect to r and 





At the leading order in slow-roll parameters, 0^ can be written as 



1 H 



3 /i 



-1/2 



(96) 



(97) 



Thus 0^ is essentially the slow-rolling parameter and it controls the strength of coupling between inflaton perturbation 
and gravitational perturbations in the bulk. We solve the coupled equations perturbatively in terms of small /3^. 



Zeroth order solutions 



At the zeroth order where = 0, the solution for u is given by 



Then Jv, becomes 



/o) = Ci(-fcr)i/V_3/2(-fcr) + C2(~fcr)i/2 J3/2(-fcr). 



V TT — fcr V TT —KT 



(98) 



(99) 



This gives the boundary condition for uj. The solution for w in the bulk subject to this condition is obtained as [I^ 



-2Ci£(-l) 

oo 

+2C2 5^(-l) 



£=0 



l^^ (sinhij'zb)Q2£(coshgz) , .-3/2 . . , n 
+ 2j sinhgzg^,(coshgz,) ("^"^ ./2.+i/2(-fcr), 

3\ sinhgz,Q2.+i(coshgz) 3^ 

2/ sinhiJzgi^_^^(coshi7zb)^ ^ ^'^^^/^^ 



where the identities 



[x) = V2^f2{-lY(2£+l)x-h2i+i/2{x), 



sin(a;) 



2^^(-l)M2£+-U-4j2,+3/2(x), 



(100) 

(101) 
(102) 



were used. 

At large scales —fcr 0, the dominant contribution comes from ^ = mode. On the other hand, on small scales 
—kr^ 00, all modes becomes comparable and we need to take into account an infinite ladder of the modes. This 
means that gravity becomes 5-dimensional at small scales. 

In practice, we must approximate the infinite sum to proceed the calculations. We first check the identity Eqs. p01|) 
and H102() to see if we can approximate the infinite summation by introducing a cut-off £c into the summation. From 
Fig. n we can see that if we increase the cut-off ic, the identity is satisfied for large —krj, i.e. on small scales. This 
implies that as long as we start from a finite time — fcr^, we can approximate the infinite ladder of the modes by 
introducing sufficiently large ic- 

Fig.Elshows the bulk solution for w{z, t) with introducing sufficiently large cut-off tc- The solution is localized near 
the brane and decays towards the horizon z ^ 00. This is a bound state that is supported by an oscillation of the 
inflaton fluctuation on the brane. This kind of bound state generally appears in coupled brane and bulk oscillators 
|23j |. A toy example is shown in Appendix. A key point here is that, in this case, the bound state is a summation of 
many different eigenstates of different eigenvalues (Eq. H10()|l '). This fact becomes crucial in the analysis of the next 
order solution. 
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FIG. 1: Dotted lines shows cos(— 2) and solid lines show the summation of Bessel functions with cut-off Ic — 20 and £c ~ 30 
respectively. 




-100 



10 -75 



FIG. 2: The zeroth order solution for uj{z,t). A brane is located at Hz = 3. 



B. First order solutions 



Now it is possible to calculate the next order equation for u^^^ 



(103) 
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where J„ describes the effect of the back reaction from the bulk perturbations. We can use the 0-th order solution to 
evaluate J„ as 



1=0 
oc 

E 

1=0 



[21+^ ^{2tH^i) (2£(2^-l)(-fcr)-ij2,+ i +2(-/cT)-ij2,+ |(-fcT)) (104) 
(-1)^ {21 + A (2£ -I- 1; Hp) (2£{2l + l)(-kT)-i J^,+ . + 2(-fcr)-^ J^f+f {-kr)) , (105) 



where 



A(n;i7/i) = - (1+ f- 



-1/2 



Q„(cosh iJzf,) 
Q^(cosh Hzi,) ' 



(106) 



The quantity A{n]Hfi) controls the amplitude of corrections to Mukhanov-Sasaki equations from the bulk over the 
change of the energy scales of the inflation. 

In order to evaluate Ju, we need to introduce a cut-off in the summation at sufficiently large £. Fig. |21 shows the 
behaviour of J„ against the change of the cut-off £c- A good feature here is that the behavior of J„ for small — fcr does 
not change even if we increase the cut-off. Thus we can reproduce the correct behaviour of Ju by a finite summation 
of modes as long as we are considering a finite time interval. 




FIG. 3: Source term Jui — z) as a function of 2 with cut-off £c ~ 20 and £c = 30 respectively. Here we take H/j, S> 1. 



1. Large scales 

On large scales — fcr ^ 0, £ = mode in Ci mode dominates, which corresponds to a = 2H^ mode. Thus 
we can approximate the infinite ladder of the modes by a single mode on super horizon scales. This indicates that, 
at large scales, gravity looks four-dimensional. Then we can easily show that J„ is suppressed for —kr and 
the Mukhanov-Sasaki equation becomes completely the same as the conventional cosmology. Thus we can show the 
conservation of the curvature perturbation TZ on large scales in the same way as conventional cosmology |2l| . 

2. Small scales 

At low energies H/fi ^ 1, A(n; H/ fi) can be approximated as 

A (n; H/fi) = il + + 1) + log(i?//i) - log 2) , (107) 
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where we assumed n is not large. Thus, the source terms is well suppressed by the term A(n; at low energies. 

However, at sufficiently small scales, large £ modes become important and the approximation 1)1071) does not hold. 
Then we could still get an effect on very sub-horizon scales {k ^ /i^^ ^ H). In this case, we need to introduce a 
large cut-off in the summation of £ and it is technically difficult to perform a calculation. 

At high energies, the amplitude of A{n;HiJ,) becomes large as Hji becomes large, but, at sufficient high energies 
Hfi oo, A{n;Hfi) becomes independent of Hfi as seen from Fig. 4. Indeed, we can obtain the asymptotic form of 
A(Z; H/ fj.) for H/fi — > oo as 

A{l;Hfx)^ i—. (108) 

n + 1 




FIG. 4: A{n;H/fj,) as a function of H/fi. 

In the following, we consider this limit. In this high energy limit, J„ is well fitted as 
2e, 



Ci(-fcr)"i/2cos(-fcT + ¥')-C'2(-fcr)-i/2sin(-/sT + (^) , (109) 



between —140 < kr < —40 where A = 0.4 and (p — 0.9. Then, the equation of motion for the first order Mukhanov 
variable is given as 

, f,.2 2 ^^ ^ (1) 1 ^^^^ ^^^^ (0) 

dr2 



<' - -7(6^ - SvK" - Mr) = 0. (110) 



By using the asymptotic behavior of the Bessel function at small scale (large —kr), the third term behaves like 
{—kr)"^ sin(— fcr), while the forth term behaves as (— fcT)"^/^ sin(— fcr). Therefore, at least at small scales, the effect 
from the bulk metric perturbations dominates the effect from the standard corrections to the de Sitter geometry. 
Thus we will neglect the third term. The general solutions are given by the linear combination of (— fcT)^/^ J3/2(— ^t) 

and (— fcr)^/2 J_3/2(— ^t). By choosing the initial conditions so that Uk{Ti) — uf\Ti) at r = r^, we find the following 
form of the solution. 



= ^(-fcT)!/^! (-fcr) + L>2(-fcr)i/2j_3(-fcT), (111) 
where Di and D2 are given by 

D2 = t\{kT'){-kT')^J^{-kT')Ju{T'). (112) 



i?i - J r d{kT'){-kT')ij_s{-kT')Mr'), 

^ Jkr, 
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For specifying the behavior of the first order Mukhanov variable, we must evaluate Di and D2. Using the asymptotic 
form for Bessel functions at small scale, Di and D2 are well approximated as 



2e . /TT 



Di ~ -^A^l- I d{kT'){-kT')-^ sm{-kT')[CiCos{~kT' + (f) - C2Sm{-kT' + (f)]. 

kTi 

D2 ~ / d{kT'){-kT')-i cos(-fcT') [Ci cos(-fcr' + (^) - C2 sin(-/sT' + ^)] . 

Then, on small scales, the first order solution is given by 

{{F{t) - Fin)) cos(-fcr) + (G(r) - G{t,)) sin(-A:T), 



(113) 



(114) 



where 



Fir) = 



G(r) 



-S 

s 
c 



f 2V^ \ (n('bEh\ 2 ^ , 
I — -jzz — 1 COS (/3 — ( G ( ^= — ) -1=^ —kT ) smiy9 



sm (/? — G 



2 



y^—kr 



COS Iy9 



2V^ 



2 ^/2V^ 

—=V—kT cos If — b ^ — I smip 

\/7r / V Jtt 



eAC2y/^ 



2V^ 



cos ip 



(115) 



(116) 



where S and G* are Fresnel functions. 

We see that the first order perturbations grows like V- 



-fcr — \J—kTi. Then if we formally take the limit —kn 00, 
the first order corrections diverge. Thus our perturbative approach breaks down. The amplitude of the zeroth order 
oscillation of inflaton fluctuations are significantly affected by the first order corrections. 

We should take care in interpreting this result for the amplitude. In a toy model of a coupled boundary and bulk 
oscillators described in Appendix A, this change of amplitude due to the first order perturbations is merely caused 
by the breakdown of the perturbative expansion. In the toy model, the coupling to the bulk oscillator just changes 
the phase of the brane oscillator. In that case we can renormalize the first-order perturbation so that the first-order 
corrections appear only in the phase of the oscillations and do not have a large effect on the amplitude. However, 
in the case of inflaton fluctuations, we cannot do this kind of renormalization. This is due to the phase ip in the 
source term of the first order equation (see Appendix A) . The phase originates from the fact that the zeroth order 
oscillation cannot be matched by a single bulk eigenmode with the same frequency as the brane oscillator and we 
need an infinite ladder of modes. Thus we can say that the effects on the amplitude from first order corrections are 
not artificial effects of our perturbative approach. 

In conventional cosmology, the amplitude of inflaton oscillations u remains constant, so we can impose initial 
conditions on any scale far inside the horizon. However, in the brane world case, the coupling to the bulk metric 
perturbations changes the amplitude of the zeroth order inflaton oscillation u, so the effect crucially depends on the 
initial conditions. In general, classically, we can also impose arbitrary initial conditions for fl. Indeed, it is always 
possible to add homogeneous solutions which satisfy the boundary condition given by 



= 0. 



(117) 



Then we find an infinite tower of massive modes starting from — 9i?^/4. Arbitrary initial conditions for SI can 
be satisfied by an appropriate summation of these massive modes. These massive modes also affect the evolution of 
inflaton fluctuations u |20l |. 

We have tried to solve the coupled equations for inflaton fluctuations and master variable directly using a numerical 
method . If we begin with the initial condition for uj given by Eq. (|100|) , the numerical solution for u well agrees with 
our perturbative solutions as long as perturbations remain valid. We have also tried using different initial conditions 
for Q and find that the effects on the amplitude of u depend on the initial conditions for lu in the bulk. 

The initial conditions for u and u must be determined by quantum theory on small scales. Thus we must quantise 
the coupled system of the inflaton fluctuations u and the master variable oj consistently. This is in contrast to the 
conventional cosmology where we can specify the vacuum for u by neglecting the gravitational effects far inside the 
horizon. This means that the assumption that the power spectrum of inflaton perturbations at horizon crossing is 
given by {H/2i:)'^ could be invalid and we may have significant effects on the amplitude of perturbations from the 
backreaction due to the bulk metric perturbations. 
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VI. CONCLUSION 

In this paper we have studied the effect of metric perturbations upon inflaton fluctuations during inflation, at first- 
order in slow-roll parameters e and r], which describe the dimensionless slope and curvature of the potential. If we 
neglect the slope and curvature of the inflaton potential then we obtain the familiar results for free field fluctuations 
in de Sitter spacetime, with a scale invariant power spectrum on large (super-horizon) scales. We take this as our 
zeroth-order result in a slow-roll expansion. 

In four-dimensional general relativity we were able to calculate corrections to the field evolution perturbatively to 
first-order in a slow-roll expansion, including linear metric perturbations. As far as we are aware this is the first time 
the slow-roll corrections have been calculated in the manner. We reproduce the familiar slow-roll corrections usually 
derived from Lyth and Stewart's exact solution to the linear perturbation equations in power-law inflation. 

On a four-dimensional brane-world, embedded in a five-dimensional bulk, there are no exact solutions for cosmolog- 
ical perturbations (for a vacuum bulk described by Einstein gravity) except for the case of an exactly de Sitter brane. 
Thus the only way to calculate slow-roll corrections is perturbatively in a slow-roll expansion. We have calculated the 
leading order bulk metric perturbations sourced by the zeroth-order inflaton fluctuations on the brane. We find that 
inflaton fluctuations support an infinite tower of discrete bulk perturbations, with negative effective mass-squared. 

Including the effect of the metric perturbations as an inhomogeneous source term in the wave equation for the 
first-order inflaton fluctuations we find that the effect of bulk metric perturbations becomes small on large scales, and 
we recover the usual result that the comoving curvature perturbation becomes constant outside the horizon. 

However at small scales (or early times for a given comoving wavelength) the effect of bulk metric perturbations 
cannot be neglected. We are able to give an approximate solution for inflaton fluctuations at high energies and on 
sub-horizon scales using a truncated tower of bulk modes. This shows that the bulk metric perturbations change the 
amplitude of inflaton field fluctuations on the brane. By including a large number of bulk modes we can model this 
effect for many oscillations, but ultimately this change of amplitude becomes a large effect leading to a breakdown of 
our perturbative analysis. 

It is not surprising in some ways that we see a large effect at small scales as these are high momentum modes which 
are expected to be strongly coupled to the bulk. Nonetheless this invalidates the usual assumption that gravitational 
effects are small far inside the cosmological horizon. It seems necessary to consistently solve for the coupled evolution 
of brane and bulk modes. We numerically tried to solve this problem and verified the validity of our perturbative 
approach as long as perturbations remain good. But it was also found that the change of the amplitude depends on 
the initial conditions for bulk metric perturbations. Detailed analysis of numerical solutions go beyond the scope of 
the present paper and they will be presented in a separate paper j23|. In order to give definite predictions for the 
amplitude of scalar perturbations in high energy inflation, we must specify the quantum vacuum state for coupled 
inflaton fluctuations and metric perturbations consistently and determine initial conditions. For this purpose, it would 
be useful to study the quantum theory of the toy model for a coupled bulk-brane oscillators in more details where we 
can consistently quantise a coupled system |23l |. 

Our result implies the possibility that the assumption that the power spectrum of inflaton perturbations at horizon 
crossing on a brane is given by (7?/27r)^ could be invalid and we may have significant effects on the amplitude of 
perturbations from the backreaction due to the bulk metric perturbations. 

Acknowledgments 

SM is grateful to the ICG, Portsmouth for their hospitality when this work was initiated and DW is grateful to the 
Maeda Lab, Waseda for their hospitality during its continuation. KK is grateful to T. Hiramatsu for the numerical 
analysis. KK is supported by the Particle Physics and Astronomy Research Council. SM is supported by the Grant- 
in- Aid for Scientific Research Fund (Young Scientists (B) 17740154). This work was also supported by PPARC grant 
PPA/V/S/2001/00544. 

APPENDIX A: TOY MODEL FOR COUPLED BULK-BRANE SYSTEM 

In this appendix, we present a simple toy model for a coupled brane and bulk oscillators. Let us consider a toy 
model for a brane field q(t) and a bulk field in Minkowski bulk, which satisfy 

4> = <P"- m'cl,, 0'(j/ = 0) = ^q. (Al) 
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We solve the equations perturbatively in terms of small /?. Without coupling, the zeroth order solution for q is given 

by 



q(0)[t) = Ci cos(Mt) + C2 sin(^t). 
If we assume m > fi, the 0-th order solution for (p is obtained as 

f3 



(Ci cos{fj.t) + C2 sin(^<)) e 



(A2) 



(A3) 



Note that the bulk field has a negative effective mass-squared and decays towards y ^ 00. This is a normalizable 
bound state supported by an oscillation of q(t) on the brane. The equation for the next order g*-^-* (t) is given by 



v(i) 



f3' 



2^r 



(Ci cos fit + C2 sin fit) . 



Including the zeroth order solution, the solution for q{t) is given by 



= Ci cos/^i + 



it sin /ii -l- C2 sin fit 



0" 



:i COS fit 



(A4) 



(A5) 



where we impose the initial condition so that q(fS) — ^'^^^(O). 

A problem is that the perturbation grows linearly in time. However, we need to be careful to interpret this growth 
of perturbations. In this toy model, we can easily find an exact solution. The corresponding exact solution becomes 



q{t) = Ci cos 



2 sm 



f3' 



Afi^/rn? 



2^r 



fi- 



-.q{t)e 



(A6) 
(A7) 



for /3 <C 1. The effect of the coupling merely changes the frequency of the brane oscillator. The origin of the linear 
instability is that the naive expansion in terms of /3 is not efhcient. Indeed, we can use 



cos(A -I- i?) = cos A cos i? — sin A sin i? ' 
sva{A + B) — sin A cos B -I- sin B cos A 



cos A — Bsin A, 
sin A + B cos A, 



(A8) 



for B 1 and expand the exact solution into Ea. ljA5|) . However, this perturbation breaks down for large t. A crucial 
difference of the inflaton fluctuations case from the toy model is that the source term for the first order equation for 
q contains a phase ip (compare Eqs. H109|l and Ulinfl to Eq. l|A4p l. Then a perturbative solution cannot be written 
into the form like Eq. (|A6|I using Eq. 1A8|I . This indicates that there could be a modification of the amplitude as well 
as the phase shift. The phase (p is originated from the fact that the brane oscillation cannot be matched by a single 
bound state (compare Eq. (|100() and Eq. (|A3|) ). Thus this is an essential difference between the toy model and the 
infiaton fluctuations case. 
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